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Radio and X-ray observations of an exceptional radio flare in the 
extreme z = 4.72 blazar GB B1428+4217 



o 
o 

(N 
X) 

IX, 

in 



> 
in 
in 

CO 
CN 
O 
VO 
O 

Oh 
6 



13 



M. A. Worsley, 1 A. C. Fabian 1 *, G. G. Pooley 2 and C. J. Chandler 3 

^Institute of Astronomy, Madingley Road, Cambridge CB3 0HA 

2 Milliard Radio Astronomy Observatory, Cavendish Laboratory, Madingley Road, Cambridge CB3 0HE 
^National Radio Astronomy Observatory, 1003 Lopezville Road, PO Box 0, Socorro, New Mexico 87801, USA 



5 February 2008 



ABSTRACT 

We report on the extreme behaviour of the high redshift blazar GB B1428+4217 at z = 4.72. 
A continued programme of radio measurements has revealed an exceptional flare in the 
lightcurve, with the 15.2 GHz flux density rising by a factor ~ 3 from ~ 140 mJy to 
~ 430 mJy in a rest-frame timescale of only ~ 4 months - much larger than any previous 
flares observed in this source. In addition to new measurements of the 1.4-43 GHz radio 
spectrum we also present the analysis and results of a target-of-opportunity X-ray observa- 
tion using XMM-Newton, made close to the peak in radio flux. Although the X-ray data do not 
show a flare in the high energy lightcurve, we are able to confirm the X-ray spectral variability 
hinted at in previous observations. GB B1428+4217 is one of several high-redshift radio-loud 
quasars that display a low energy break in the X-ray spectrum, probably due to the presence 
of excess absorption in the source. X-ray spectral analysis of the latest XMM-Newton data is 
shown to be consistent with the warm absorption scenario which we have hypothesized pre- 
viously. Warm absorption is also consistent with the observed X-ray spectral variability of the 
source, in which the spectral changes can be successfully accounted-for with a fixed column 
density of material in which the ionization state is correlated with hardness of the underlying 
power-law emission. 
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1 INTRODUCTION 

Blazars are some of the most extreme extragalactic sources known: 
powerful, radio-loud active galactic nuclei (AGN) in which the 
relativistic jet is aligned with the observer's line of sight. Their 
spectral energy distribution (SED) is dominated by beamed, non- 
thermal continuum emission from the jet and has a character- 
istic double-humped shape, with synchrotron radiation from the 
jet dominating at low energies, and Compton up-scattered emis- 
sion dominating at high energies. In BL Lac objects (the lower- 
luminosity class of blazars), the beamed synchrotron emission can 
peak in the soft X-ray regime, with the second peak dominating 
the TeV range. For flat-spectrum radio-quasars (FSRQs; the most 
luminous blazars), the synchrotron radiation tends to peak in the in- 
frared, with the hi gh energy emission peaking in the 7-ray regime 
jFossatietaflll998h and hard (T ~ 1.5) power-law in the X-ray 
range. 

In addition to their extreme luminosity and characteristic SED, 
blazars often show spectacular broad-band variability, in both ab- 
solute brightness and spectral shape. Flares in flux of several tens 
of per cent have been observed in the radio to TeV regimes in sev- 
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eral sources, on timescales of minutes to years depending on wave- 
length. Blazars show strong, core-dominated and highly-polarised 
radio emission, strongly supportive of the relativistic jet interpreta- 
tion. Blazars are usually unresolved although very long baseline 
interferometry (VLBI) can often reveal superluminal motion on 
milliarcsecond scales within the nucleus, again consistent with the 
presence of a relativistic jet. 

GBB1428+4217, at a redshift of z = 4.72, shows an 
unexpected flattening of the X-ray power-law continuum at 
soft X-ray energies. This was first seen in ROSAT a n d Bep- 
poSAX observations iBoller et alj l2000t IFabian et al] k oOla) 
and su bsequently confirmed with XMM-Newton Iworsiev^t^l] 
12004 ah . A flattening in the X-ray spectral index below ob- 
served energies of ~ 1 keV has been observed in several other 
high-redshi ft, radio-loud sou rces, e.g.: RBS315 at z = 2.69 
{Pico ncelli & Guainazzil l2005h: PMN J0525 — 3343 at 2 = 4.41 
iFabian et alj bOOlbl IWorslev et alj 12004th : RX J1028.6-0844 
at z = 4.2 8 lYuanet alJ l200(j l2005d) : PMN J1451 — 1512 at 
z = 4.76 lYuari et all l2005bh : PKS B 1830-21 1 at 2 = 2.51 
| de Rosa et all l2005h: and PKS 2126-0158 at z = 3.27 
iFiore et alj l2003t iFerrero & Brinkmannl [2003). Although low 
energy miscalibrations i n ASC A were initially proposed as an 
explanation iGrupe et aljl2004h . the large number of high-quality 
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XMM-Newton observations have shown beyond doubt that the ef- 
fect is real in several objects; that said, a number of XMM-Newton 
observations have also shown that several other high-redshift 
and radio-loud quasars show no evidence of sp ectral flattening 
iFiore et all2003l;lGrupe et alJl2004EoOrilPiconcelli et all2005h . 
A trend toward spectral flattening with increasing redshift has also 
been s een in several studi e s of radio-loud obi ects (e . gJCagrji^yilJ 
19971 IFiore et alJ Il998t beeves & Turned 1200(1 iBassett et ail 



2004) 



The origin of the spectral depression is not fully under- 
stood, and, whilst the spectral shape of GBB 1428+4217 (and 
other sources) has been repeatedly characterised using an absorbed 
power-law model, there is not yet sufficient evidence to be certain 
that intrinsic absorption is the solution. Both Galactic and inter- 
galactic absorption can be effectively discarded as viable solutions 
given the high column densities required and, critically, by the fact 
that the phenomenon h as not been observed i n radio-quiet objects 
at comparable redshifts fchemmer et alfeOOllvignali et all2 005). 
which would be expected to be similarly affected by absorption 
systems lying along the line of sig ht I O'Flahertv & Ja kobserll997t 
Fabian et al. 2001b; Worslev et al. 2004b). The remaining possibil- 
ities are absorption intrinsic to the source or an underlying spectral 
break in the power-law continuum; however, neither explanation 
can yet be ruled-out in any of the cases of soft X-ray flattening. 

As discussed bv lFabian et aljl2001bl) . an underlying break in 
the continuum could arise in a number of situations. The hard X- 
ray component of the blazar SED is believed to be due to inverse 
Compton scattering of seed photons by electrons in the relativistic 
jet. A spectral break can arise through a low-energy cut-off in the 
electron population or by a sharply peaked seed photon distribution. 
Both mechanisms are rather speculative given our limited knowl- 
edge of the processes which are important in blazar emission and 
difficult to reconcile with the sharpness of the X-ray break, which 
occurs over a range of only a few keV in the rest-frame. An in- 
trinsic absorber, with a column density of TVh ~ 10 22 -10 23 cm -2 , 
remains the most self-consistent explanation for the soft X-ray flat- 
tening although very little is known about the nature or origin 
of the material, which, given the lack of the effect in radio-quiet 
sources, seems to be closely linked to the presence of jets in these 
young AGN. X-ray data alone have been unsuccessful in break- 
ing the degeneracy between a cold (neutral) absorber or one that 
is warm (ionized). In some cases, e.g. for PMN J0525 — 3343 and 
RX J1028.6— 0844, optical observations have been able to rule- 
out cold absorption from the lack of a Lyman-limit system at the 
redshift of the quasar by the presence of emission at rest-frame 
wavelengths shortward of 912 A. The lack of such a system in 
an object means the column density in neutral hydrogen must be 
N-g.1 < 3 x 10 17 cm -2 , significantly less than the X-ray implied 
JVhi ~ 10 22 ~10 23 cm -2 . 

In this paper we present the analysis and results from con- 
tinued radio monitoring of GB B 1428+4217, in particular, the de- 
tails of an exceptional flare in the lightcurve which has seen the 
15.2 GHz flux increase by a factor of ~ 3 in only ~ 4 months in 
the rest-frame of the source. We also report the findings of a recent 
XMM-Newton target-of-opportunity observation performed close to 
the peak in the radio lightcurve and discuss the X-ray luminosity 
and spectral evolution displayed by this source. A broad-band radio 
spectrum is also presented. We ado pt Hp = 7 1 km s~ Mpc - , 
Q M = 0.27 and Q A = 0.73 ISpergel et all2003l) . Errors are quoted 
at the lcr (68 per cent) level unless stated otherwise. Also note that 
any X-ray fluxes quoted are always corrected for Galactic absorp- 



tion but are not corrected for any intrinsic absorption unless stated 
otherwise. 



2 PREVIOUS OBSERVATIONS 

GBB1428+4217 was discovered bv lHook & McMahorJ fl993) as 
part of survey to identify FSRQs with red optical counterparts 
(GBB1428+4217 has R = 20.9 with B-R> 3.40). It is at a 
redshift of z — 4.72, making it one of the most distant X-ray 
sources known. The blazar has an inverted radio-spectrum, with a 
1.4-9 GHz spectral index of ai.4_9 — — 0.4(/„ oc v~ a ). Imaging 
revealed a core-dominated sourc e with a hig h level of polarisation 
(5.5 per cent at 1.5 GHz; Con don et alJll998l) . The optical spec- 
trum reveals broad emission lines due to Lya, N V, the Si IV/O IV] 
blend and C IV iHook & McMahorll998l) . We will first summarise 
the findings of previous X-ray and radio observations, before we 
move on to present the latest results and a joint discussion of the 
spectrum and variability. 

2.1 X-ray spectrum 

An X-ray counterpart was found in an archival ROSAT Posi- 
tion Sensitive Proportional Counter (PSPC) observation taken 
in 1992. Subsequent observations we re performed wi th the 
ROSAT High Resolution Imager (HRI; iFabian et alJ \l99% and 
ASCA I Fabian et al. 1998). The inferred isotropic X-ray lumi- 
nosity exceeded 10 47 erg s _1 which suggested highly beamed 
emission. The 0.5-10 keV data were consistent with a hard 
power-law spectrum, with a photon index of T = 1.29 ± 0.05, 
and an upper limit to any intrinsic absorption at the red- 
shift of the quasar of 4.5 x 10 22 cm -2 (90 per cent confi- 
dence); however, the authors noted that the SIS instrument alone 
did point to some excess absorption - a column density of 
(1.3 ± 0.6) x 10 23 cm" 2 forT = 1.43 ± 0.12. Excess absorption 
was confirmed with a further ROSAT PSPC observation, which 
found a power-law index of T = 1.4 ± 0.2 and a column den- 
sity of ATg^_=jUj2±0 ; 28) x 10 22 cm -2 (90 per cent errors) at 
z — 4.72 jBoller et al J200d) . Further observations with BeppoSAX 
also measured Y — 1.45 ± 0.10 power-law and exce ss abso rption 
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Warm absorption was first suggested by Fab ian et alJ 1200 lal) 
following the BeppoSAX observations. In addition to mod- 
elling the results of that observation, the authors also fit- 
ted the previous ASCA and ROSAT data. A marginally con- 
sistent set of parameters were found with a column density 
of (2-3) x 10 23 cm -2 and an ionization param eter £ ~ 300. 1 
XMM-Newton observations I Worslev et ai]l2004iJ) confirmed the 
presence of soft X-ray spectral flattening, with a column 
density of A^ old ~ (1.4-1.6) x 10 22 cm" 2 if cold, or up to 
JV™ m ~ 10 23 cm -2 with an ionization parameter £ ~ 100, if due 
to a warm absorber. Interestingly, the XMM-Newton data showed 
a distinct change in the spectral index of the source, both from 
the earlier observations as well in the ~ 7 d (rest-frame) between 
the two XMM-Newton exposures. The first observation indicated 



1 IFabian et alj 12001 al) quote Ci-1000 Rydi denned as 
^1-1000 Ryd/nuR 2 , where L 1 _ 1000 Ryd is the 1-1000 Rydberg 
luminosity, njj is the space dens ity of H I and R is t he distance from the 
source. For GB B1428+4217, the lFabian etai]<2001al) definition is a factor 
~ 2 higher than the £ = L2— 10 keV /ih^ 2 definition used in this paper. 
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r ~ 1.9, whilst the second found V ~ 1.75. Whilst the source had 
long been known to show extreme luminosity variations, this was 
the first definitive example of spectral changes. 

2.2 X-ray variability 

No clear variability was noted between the original ROSAT 
PSPC detection in 1992 and the ROSAT HRI and ASCA ob- 
servations during 1996 and 1997: all were consistent with a 
0.1-2.4 keV soft X-ray flux 2 of ~ 10" 12 e rg cm" 2 s" 1 , with 
no significant evidence for any differences I Fabia n et alJll997L 
1998). However, a series of four ROSAT HRI observations from 
1997 December 12 until 1998 January 23 showed substantial 
X-ray variability, with the 0.1-2.4 keV flux increasing from 
(0.6 ± 0.1) x 10" 12 to (1.1 ± 0.1) x 10" 12 erg cm" 2 s" 1 be- 
tween the third and fourth exp osure - a period o f only 2.4 d 
in the rest-frame of the source <Fabianetallll999h . Subsequent 
ROSAT PSPC observations then identifie d 25 per cent va riabil- 
ity on a rest-frame timescale of ~ 1.8 h IB oiler et all 2000). Fur- 
thermore, the mean 0.1-2.4 keV flux during the observation was 
(2.7 ± 0.6) x 10" 12 erg cm" 2 s" 1 , 3 times the previous ROSAT 
and ASCA measurements in this band. A cautionary note must 
be attached to these PSPC measurements, which were amongst 
the final observations of the ROSAT satellite which experienced 
some detector anomalies and the development of a large gain 
hole in the PSPC instrument, although careful data reduction tech- 
niques (refer to lBoller et ai]|200d) should mean the results were 
not sig nificantly affected. The BeppoSAX-menswsd l Fabi an et alj 
[2001bb 2-10 keV flux and spectral shape imply a 0.1-2.4 keV 
flux of ~ 10 -12 erg cm -2 s" 1 , indicating a decrease back to 
the level seen prior to the ROSAT PSPC observation. The two 
XMM-Newton exposures of 2002 December 09 and 2003 Jan- 
uary 17 measured 0.1-2 keV fluxes of (1.17 ±0.04) x 10" 12 
and (0.85 ± 0.02) x 10" 12 erg cm" 2 s" 1 respectively, consistent 
with other measurements but revealing a ~ 30 per cent decrease in 
the 7 d (rest-frame) between the XMM-Newton observations. 

In the hard X-ray band, ASCA and BeppoSAX l Fabi an et alJ 
Il998ll2001bl) measured consistent 2-10 keV fluxes of 2.5 x 10" 12 
and 2.76 x 10" 12 erg cm" 2 s" 1 . The two XMM-Newton expo- 
sures measured 2-10 keV fluxes of (1.21 ± 0.06) x 10" 12 and 
(1.08 ± 0.03) x 10~ 12 erg cm -2 s _1 respectively. These were 
~ 2-2.5 times lower than the earlier measurements and indicating 
a 10 per cent decrease in only 7 d (rest-frame) between the XMM- 
Newton observations (note that this decrease is much less than the 
~ 30 per cent fall seen for the soft X-ray band because of the dis- 
tinct change in the spectral shape of the source between the two 
observations: see Section l2~TV We will present soft and hard band 
X-ray lightcurves shortly, following the details and results of the 
new XMM-Newton measurements. 

2.3 Radio properties 

Fabia n~ alJ 1 1997b collated various radio measurements for 
GBB1428+4217 and present the SED of the source, noting the 
flat radio spectrum, high level of polarisation and indications of 
5 GHz variability. An intensive monitoring programme was car- 
ried out between 1998 April and December using the Ryle Tele- 
scope at Cambridge at 15.2 GHz. The programme formed part of 

2 To reiterate, all X-ray fluxes have been corrected for Galactic absorption 
only, and not for any apparently intrinsic absorption. 



Table 1. VLA measurements of the GB B 1428+4217 radio spectrum from 
1.425 to 43.4 GHz as measured on 2005 October 14 (MJD 53657), close 
to the peak in the radio lightcurve. The 15.2 GHz data are from the Ryle 
Telescope on 2005 October 12. Ryle Telescope flux densities are of Stokes' 
I + Q but the VLA measurements are of Stokes' I. The flux densities of 
the calibration source 3C286 are also given. 



Frequency 
(GHz) 


GBB1428+4217 
flux density (mJy) 


3C286 calibrator 
flux density (Jy) 


1.425 


171 ± 9 


14.70 


4.86 


273 ± 14 


7.47 


8.46 


348 ± 17 


5.20 


15.2 


380 ± 38 


3.50 


22.46 


367 ± 37 


2.52 


43.34 


290 ± 29 


1.45 



a radio and X-ray variability study, although the radio observations 
were somewhat later than the ROSAT HRI exposures carried out in 
1997 December and 1998 January jFabian et alii 19991) . The Ryle 
data revealed significant 15.2 GHz variability, with a flux increase 
from ~ 100 mJy to ~ 140 mJy in a rest-frame period of only 
~ 20 d. The flare ended abruptly mid-way through the observing 
programme, with the flux decreasing back to ~ 105 mJy within a 
fraction of a day in the rest-frame. The lightcurve also shows strong 
variability at the 15 per cent level with a rest-frame timescale of 
only ~ 2 d. 



3 NEW RADIO OBSERVATIONS 

The Ryle telescope has been used to measure the 15.2 GHz ra- 
dio flux of GBB 1428+4217 on a regular basis since the monitor- 
ing programme in 1997-8 (for details of the observing method see 
IPoolev & FendeJ 1997b . The lightcurve is shown in Fig.0a). Mea- 
surements after the 1997-8 programme, although sparse, indicated 
no significant behaviour beyond the variability that had already 
been observed until that point. The flux decreased to only 70 mJy 
before starting to climb. In less than 6 months (rest-frame), the flux 
peaked at ~ 210 mJy before dropping suddenly to 140 mJy. The 
source then began to steadily grow in flux - the start of a flare of 
a enormous magnitude that has continued to the present. By 2005 
August, the 15.2 GHz flux had reached ~ 430 mJy, a factor ~ 3 
increase in only ~ 4 months (rest-frame). The lightcurve continues 
to show intra-day variations of ~ 10 per cent, with several 'micro- 
flares' visible in the primary flare. 

A radio spectrum was obtained on 2005 October 14 with 
the Very Large Array (VLA) at 1.425, 4.86, 8.46, 22.46 and 
43.34 GHz, and almost simultaneously at 15.2 GHz with the Ryle 
Telescope, shortly after the peak flux density at 15.2 GHz (flux 
densities at this frequency are of Stokes' I + Q\ all VLA measure- 
ments are of Stokes' J). The data and spectrum are given in TableQ 
and Fig. [2] respectively. The absolute flux calibration for the VLA 
data was done using 3C286 and the measured flux densities at each 
frequency are given in Table Q The flux measurements include an 
uncertainty in the calibration of 5 per cent for the three lowest fre- 
quencies and 10 per cent for the remainder. 
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Figure 1. (a) The 15.2 GHz (rest-frame 86 GHz radio lightcurve of 
GBB1428+4217 as measured by the Ryle Telescope in Cambridge from 
1998 March 20 (MJD 50923) until 2005 November 01 (MJD 53623). (b) 
The soft 0.1-2.4 keV (rest-frame ~ 0.6-14 keV) X-ray lightcurve for the 
same time period, (c) The hard 2-10 keV (rest-frame ~ 11-57 keV) X-ray 
lightcurve. X-ray Measurements and errors have been taken directly from 
the reference paper or, where possible, estimated from the spectral parame- 
ters and other details. References (in order of observation date): (a) ROSAT 
HRI IFabian et aljli 



19971) : (b) ASCA IFabian et alj|l998l): (c) ROSAT HRI 



Fabian et al. 1999); (d) ROSAT PSPC i Boiler et al. 2000); (e) BeppoSAX 
Fabian et all2001al) : (f ) XMM-Newton rev. 5 49 IWorslev et all2004al) : (g) 

XMM-Newton rev. 569 IWorslev et all2004al) : (h) XMM-Newton rev. 1005 

(this paper). 



4 NEW XMM-NEWTON OBSERVATIONS 
4.1 Data reduction 

Following the extraordinary increase in 15.2 GHz flux, XMM- 
Newton was used to perform a target-of-opportunity (ToO) X-ray 
observation co-incident with the radio brightening. Various details 
of the observation, which was carried out during XMM-Newton rev- 
olution 1005, are given in Table [2] The European Photon Imaging 
Camera (EPIC) pn, MOS-1 and MOS-2 instruments were active. 
Data were reduced using the XMM-Newton Science Analysis Sys- 
tem (SAS) version 6.4.0-6.5.0. Periods of background flaring were 
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Figure 2. A vF v plot showing the radio spectrum of the blazar using data 
from the VLA and Ryle Telescope. The data are given in TableQ 



excluded using a lightcurve extracted from the whole pn/MOS chip 
in each case. Raw events were filtered in the recommended 3 way. 

Spectra were produced from events extracted from circular 
regions centred on the source. Radii of approximately 45 and 
75 arcsec were used for the pn and MOS, respectively. Background 
spectra were formed using events from nearby source-free regions 
on the same chip as the source. The regions used had respective 
radii of approximately 130 and 170 arcsec for the pn and MOS. 
The SAS tasks ARFGEN and RMFGEN were used to generate the 
appropriate response matrices. Background-corrected spectra were 
produced for each instrument and grouped to a minimum of 20 
counts in each energy channel. 



4.2 X-ray lightcurve 

Figs. 0b) and 0c) show X-ray lightcurves to match the ra- 
dio data shown in Fig. Qa). Measurements of the X-ray flux 
in the soft 0.1-2.4 keV and hard 2-10 keV bands (corre- 
sponding to rest-frame energies of ~ 0.6-14 and ~ 11-57 keV, 
respectively) are shown. The measurements are taken from 
the various observations described in Section 12.21 and the 
new XMM-Newton data (2005 June 05). This observation had 
0.1-2.4 and 2-10 keV fluxes of (0.70 ± 0.02) x 10~ 12 and 
(1.27 ± 0.04) x 10~ 12 erg cm' 2 s _1 , respectively. The X-ray 
lightcurve illustrates the strong variations in flux that have been 
seen throughout the various observations of the source (see Sec- 
tion |2]2| for details and references). The latest XMM-Newton data 
do not show any evidence of variability within the observation; 
however, the 0.1-2.4 and 2-10 keV fluxes are different to those 
measured in the previous XMM-Newton observations. 

There appears to be no obvious correlation between the X- 
ray and radio lightcurves, although the sparse sampling of the 
lightcurves does not allow us to rule anything out. Interestingly, 
one would not necessarily expect a correlation with zero time-lag, 
particularly since blazar emission is dominated by beamed radia- 
tion from a relativistic jet. Continued monitoring in both radio and 



A (FLAG==0) && (PATTERN<=4) && (PI IN [150, 15000]) and 
(&XMMEA_EM) && (PATTERN< = 12 ) && (PI IN [200, 10000]) for 
the pn and MOS cameras, respectively. See the XMM-Newton Science 
Operations Centre document XMM-SOC-CAL-TN-0018. 
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Table 2. Details of the ToO XMM-Newton observation from revolution 1005 along with the details of the previous exposures of revolutions 549 and 569. 



Revolution Date MJD Mode Filter Good Exposure Times (ks) Good counts 













pn 


MOS-1 


MOS-2 








549 


2002 Dec 09 


52617.3 


Full frame 


Thin 


2.6 


4.1 


4.3 


2162 


996 


979 


569 


2003 Jan 17 


52656.7 


Full frame 


Thin 


11.5 


14.2 


14.2 


6797 


2577 


2541 


1005 


2005 Jun 05 


53526.3 


Full frame 


Thin 


12.9 


15.5 


15.9 


7222 


2604 


2657 



X-ray wavebands over the coming years may reveal correlations 
which are not clear with the present data. 

4.3 X-ray spectrum 

XSPEC version 11.3 was used for spectral fitting, which was done 
for the pn, MOS-1 and MOS-2 data jointly. A simple power-law 
model 4 , fitted at energies > 1 keV, is a good match to the spec- 
trum, with a reduced chi-squared of xt = 1-024 for v = 316 de- 
grees of freedom (the probability of the model not representing the 
data, 1 — Pr(Ho), is ~ 63 per cent and there is therefore no sig- 
nificant evidence to suggest the model is unacceptable). However, 
when the same model is applied to the full energy range, the fit 
is unacceptable (at the 99.998 per cent level), with xlss ~ 1-359. 
Re-fitting the simple power-law model to the full energy range 
also leads to an unacceptable fit (at the 99.7 per cent level), with 
xl ss = 1-187. 

Soft X-ray spectral flattening is once again evident, and we use 
a powe r-law plus intrinsic co ld absorption model (using ZWABS in 
XSPEC; lMorrison & McCammonlll983l) . A joint analysis was per- 
formed with the earlier XMM-Newton data taken in 2002 Decem- 
ber and 2003 January (XMM-Newton revolutions 549 and 569), re- 
spectively; the details are included in Table |2] (for further infor- 
mation on these data refer to lWorslev et all2004al) . Table[3]shows 
the various results. Allowing the power-law photon index F, the 
ZWABS column density A^H° ld , and the normalization to vary be- 
tween observations can yield good fits to the revolution 549 and 
569 data, although the 1005 data is less satisfactory. Fixing F and 
Nft ld between the observations results in an unacceptable fit (at 
the 99.994 per cent level), with X1121 = 1-170. Alternatively, al- 
lowing the A^H 0ld to vary whilst T is fixed is also unacceptable, 
with X1119 = 1-111 (at the 99.5 per cent level): there is clear spec- 
tral variability in the source (a discussion of which we will come to 
shortly). Interestingly, fixing jV^ arm but allowing T to vary leads to 
a statistically acceptable fit, with X1119 = 0.999. All three data-sets 
are explicable with variations in the underlying power-law slope, 
with a fixed column density of (1.4 ± 0.1) x 10 22 cm~ 2 of intrin- 
sic cold absorption. This model is not a good fit to the revolution 
1005 data on its own (x1s4 ~ 1-095; unacceptable at the 99.7 per 
cent level), something which may suggest that the spectrum is more 
complex, perhaps involving warm rather than cold absorption. 

Warm absorption is also a very plausible candidate for the 
spectral shape of the source. Table |4| shows the results of various 
power-law plus warm absorption mod el fits to the data. The pho- 
toionization code CLOUDY I Ferland et al. 1998) was used to gen- 
erate the warm absorber models. The ionization parameter £ is de- 
fined from £ = L/nuR 2 ', where L is the rest-frame 2-10 keV lu- 
minosity, nu is the absorber space density and 7? is the distance 

4 Ga lactic cold absorption of JV^ al = 1.40 X 10 20 cm" 2 lElvis et alj 
1 1994 . applied using WABS in XSPEC iMorrison & McCammorill983ft r7s 
included in all spectral models 



from the source. As with the cold absorption models we tested var- 
ious combinations of free and fixed parameters in the analysis. Al- 
lowing r to vary along with both N^ ld and the ionization param- 
eter £ (this is a fixed £ = in the cold absorber models) results 
in a good fit with X1114 — 0.985. Fixing T results in a poorer fit 
with X1127 — 1.059, although the evidence for an unacceptable fit 
is only at the 91.7 per cent level. Allowing T and £ to vary whilst 
maintaining a fixed TV™ 1 " 1 provides a very good description of the 
data, with X1126 = 0.987 and reasonable values for the fits to the 
individual data-sets. A fixed column density of absorption with ion- 
isation state and underlying power-law variations seems the more 
physically plausible situation and is also consistent with the X-ray 
spectral variation which we will describe next. 



5 X-RAY SPECTRAL VARIABILITY 

The X-ray lightcurves (Figs. Qa) and 0b)) show strong ampli- 
tude variations but also show differences between the 0.1-2.4 and 
2-10 keV behaviour. The (0.1-2.4)/(2-10) keV flux ratio was 
~ 0.5 and then ~ 0.4 during the ASCA and BeppoSAX observa- 
tions. The ratio was much larger during the later XMM-Newton ob- 
servations, with more evidence of spectral variation between the 
revolution 549 and 569 observations, with the ratio falling from 
0.97 ± 0.06 to 0.78 ± 0.03. The revolution 1005 data indicate a 
value of 0.55 ± 0.02, close to that seen in the earlier observations. 
Here, we have used a simple soft/hard band flux ratio to charac- 
terise the spectral variation but the evolution can also be well- 
described by a change in the underlying power-law slope with a 
fixed column density of absorption. Fig. [3] shows the X-ray spec- 
tra of the three XMM-Newton observations and clearly indicates a 
hardening of the power-law photon index in moving from the revo- 
lution 549 to the 569 and 1005 data. 

Fig. shows the confidence contours in the power-law slope 
r and the column density of an intrinsic neutral absorber, 

as measured by the ROSAT PSPC, BeppoSAXand the three XMM- 
Newton observations. Interestingly, the fixed- N^ ld but varying-r 
scenario (Table |5Jd)), is marginally consistent with all the data- 
sets. The spectral variations can be explained through variations 
in the photon index with a fixed absorber column density of 
~ 1.4 x 10 22 cm -2 . These variations, however, do not seem to 
be correlated with changes in the normalization of the power-law 
(which is presumably indicated by the 2-10 keV lightcurve in 
Fig.EJc)). 

The spectral variation can also be explained successfully in 
the warm absorber scenario. The most self-consistent situation is 
that of a fixed-A^H arm but varying-r and varying-^, i.e. the col- 
umn density of the absorber remains fixed whilst the ionization 
state increases with the hardness of the irradiating power-law; 
indeed, £ seems to increase with the hardening T from revolu- 
tions 549, 569 and 1005 (Table [4}- This trend can also be visu- 
alised in Fig. [5] which shows the X-ray spectrum hardening be- 
tween the data from XMM-Newton revolutions 549, 569 and 1005. 



6 M. A. Worsley et al. 



Table 3. Spectral fits to apower-law plus intrinsic cold absorption model (inclusive of Galactic cold absorption of N^ 3,1 = 1.40 X 10 20 cm -2 ). The fits use 
different combinations of parameters which were fixed to be the same for the separate observations, and those which were free to vary independently. The xl 
and degrees of freedom v is shown for the total fit and also when the model is applied to the individual data set alone. 



Revolution Y N™ ld Total Individual 

(10 22 cm" 2 ) xl v xl v 

(a) T and N^ M independent: 

549 1.89 ±0.04 1.61 ± 0.28 0.831 185 

569 1.76 ±0.02 1.67 ±0.16 0.991 1117 0.946 454 

1005 1.49 ±0.02 0.95 ±0.15 1.081 484 

(b) T and A^ old fixed: 

549 * 1.189 185 

569 1.67 ±0.01 1.40 ±0.10 1.170 1121 1.005 454 

1005 1.308 484 

(c) T fixed; 7V™ ld independent: 

549 0.51 ±0.17 1.031 185 

569 1.67 ±0.01 1.15 ±0.13 1.111 1119 0.989 454 

1005 2.12 ±0.15 1.245 484 

(d) T independent; N™ ld fixed: 

549 1.86 ±0.27 0.835 185 

569 1.73 ±0.02 1.4 ±0.1 0.999 1119 0.953 454 

1005 1.53 ±0.02 1.095 484 



Table 4. Spectral fits to a power-law plus intrinsic warm absorption model (inclusive of Galactic cold absorption of 7V^ al = 1.40 X 10 20 cm -2 ). The fits 
use different combinations of parameters which were fixed to be the same for the separate observations, and those which were free to vary independently. The 
X 2 and degrees of freedom v is shown for the total fit and also when the model is applied to the individual data set alone. 



Revolution Y log 10 (A r ^ arm /cm -2 ) log 10 (£/erg cm s _1 ) Total Individual 

xl v xl v 

(a) T, AT™™ md £ independent: 

549 1.89 ±0.06 22.34 ± 0.21 < -1 0.834 185 

569 1.77 ±0.04 22.37 ±0.15 -0.04 ± 1.04 0.985 1114 0.945 454 

1005 1.55 ±0.04 22.90 ±0.67 2.03 ± 0.44 1.061 484 



(b) T fixed; 7V™ arm and £ independent: 

549 ' 21.82 ±0.50 < -1 0.992 185 

569 1.70 ±0.03 22.24 ±0.13 -0.04 ± 0.96 1.059 1127 0.969 454 

1005 > 23 1.84 ±0.26 1.153 484 



(c) T and § independent; 7V™ arm fixed: 

549 1.91 ±0.06 -0.96 ±0.30 0.836 185 

569 1.77 ±0.02 22.34 ± 0.04 -0.18 ± 1.04 0.987 1126 0.946 454 

1005 1.51 ±0.02 1.44 ±0.15 1.071 484 



A fixed log 10 (A r H arm /cm -2 ) ~ 22.5 absorber is consistent with 
all the data-sets, with the underlying power-law softening from 
F ~ 1.4 to ~ 1.9, and then hardening back to ~ 1.5 during 
the BeppoSAX and the three XMM-Newton observations. A trend 
is also seen in the ionization parameter, which decreases from 
log 10 (£/erg cm s -1 ) ~ 1.5 to ~ 0.8, and then increases to 1.6. 



6 CONCLUSIONS 

• We have observed an exceptional flare in the 15.2 GHz radio 
lightcurve of GB B 1428+4217. The flux density has grown by a 
factor of ~ 3 in a rest-frame timescale of only ~ 4 months. The 
radio lightcurve also shows intra-day variations on the ~ 10 per 
cent level, consistent with previous observations. 

• Target-of-opportunity X-ray observations with XMM-Newton 



do not see any evidence of a flare in the X-ray lightcurve, although 
the source continues to show ~ 10-25 per cent variability on rest- 
frame timescales of hours to days. 

• The ToO observations are able to confirm the unusual 
X-ray spectrum of the source, which is well-described 
by a hard (r ~ 1.4-1.9) power-law with a large amount 
(~ 10 22 -10 23 cm" 2 ) of intrinsic absorption (refer to 
Worsl ev et all 12004 sk . The present data are unable to dis- 
criminate between a cold (neutral) or warm (ionized) absorption (a 
third possibility, of a spectral break in the power-law continuum, is 
discussed in detail in our previous paper). 

• Previous XMM-Newton observations hinted at spectral varia- 
tions in the source. These have been confirmed with the ToO data, 
which reveals that the photon index of the X-ray power-law has 
hardened from T ~ 1.7-1.9, to F ~ 1.5. Interestingly, there is no 
clear evidence for a change in the column density of the absorption 



Radio and X-ray observations of an exceptional radio flare in the extreme z = 4.72 blazar GB B1428+4217 
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Figure 3. Unfolded X-ray spectra, in vF v space, from the XMM-Newton 
revolution 549 (open black circles), 569 (solid grey diamonds) and 1005 
(solid black stars) observations. A power-law plus intrinsic cold absorption 
model was used to unfold the spectra and is shown as a solid line in each 
case. 
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Figure 5. Confidence contours (68, 90 and 99 per cent) in the photon in- 
dex r and the warm absorber column density A r ^ arm . Contours are also 
shown for the BeppoSAX data I Fabian et al. 2001 ^) and the da ta from XMM- 
Newton revolutions 549, 569 and 1005 jWorslev et all2004al this paper). A 
fixed- jV^ arm warm absorber can account for all four observations through 
strong variations in the power-law slope and ionization parameter. 
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Figure 4. Confidence contours (68, 90 and 99 per cent) in the photon index 
T and the c old absorber column density JV£ old . The BeppoSAX and ROSAT 
PSPC data <Fabian et aj]|2001atlBoller et alj|200fj) are sho wn, along with 
results from XMM-Newton revolutions 549, 569 and 1005 i Worslev et al. 
l2004it this paper). The data are all marginally consistent with a constant 
N™ M model in which V varies widely, although the column densities im- 
plied by the revolution 569 and 1005 data disagree at the 90 per cent level. 



in either the cold or ionized absorber scenario. The spectral vari- 
ations, including all previous data from other telescopes, are con- 
sistent with a scenario in which the column density of an ionized 
absorber remains fixed whilst the ionization state increases with the 
hardening of the power-law emission. 
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